Fundamentals  of  Signal  Processing  for  Phased  Array  Radar 


Dr.  Ulrich  Nickel 

Research  Institute  for  High-Frequency  Physics  and  Radar  Techniques  (FHR), 
Research  Establishment  for  Applied  Science  (FGAN), 

53343  Wachtberg,  Germany 

nickel@fgan.de 


ABSTRACT 

This  section  gives  a  short  survey  of  the  principles  and  the  terminology  of  phased  array  radar. 
Beamforming,  radar  detection  and  parameter  estimation  are  described.  The  concept  of  subarrays  and 
monopulse  estimation  with  arbitrary  subarrays  is  developed.  As  a  preparation  to  adaptive  beam  forming, 
which  is  treated  in  several  other  sections,  the  topic  of  pattern  shaping  by  deterministic  weighting  is 
presented  in  more  detail. 


1.0  INTRODUCTION 

Arrays  are  today  used  for  many  applications  and  the  view  and  terminology  is  quite  different.  We  give  here 
an  introduction  to  the  specific  features  of  radar  phased  array  antennas  and  the  associated  signal  processing. 
First  the  radar  principle  and  the  terminology  is  explained.  Beamforming  with  a  large  number  of  array 
elements  is  the  typical  radar  feature  and  the  problems  with  such  antennas  are  in  other  applications  not 
known.  We  discuss  therefore  the  special  problems  of  fully  filled  arrays,  large  apertures  and  bandwidth.  To 
reduce  cost  and  space  the  antenna  outputs  are  usually  summed  up  into  subarrays.  Digital  processing  is 
done  only  with  the  subarray  outputs.  The  problems  of  such  partial  analogue  and  digital  beamforming,  in 
particular  the  grating  problems  are  discussed.  This  topic  will  be  reconsidered  for  adaptive  beamforming, 
space-time  adaptive  processing  (STAP),  and  SAR. 

Radar  detection,  range  and  direction  estimation  is  derived  from  statistical  hypotheses  testing  and 
parameter  estimation  theory.  The  main  application  of  this  theory  is  the  derivation  of  adaptive 
beamforming  to  be  considered  in  the  following  lectures.  In  this  lecture  we  present  as  an  application  the 
derivation  of  the  monopulse  estimator  which  is  in  the  following  lectures  extended  to  monopulse  estimators 
for  adaptive  arrays  or  STAP. 

As  beamforming  plays  a  central  role  in  phased  arrays  and  as  a  preparation  to  all  kinds  of  adaptive 
beamforming,  a  detailed  presentation  of  deterministic  antenna  pattern  shaping  and  the  associated  channel 
accuracy  requirements  is  given. 


2.0  FUNDAMENTALS  OF  RADAR  AND  ARRAYS 

2.1  Nomenclature 

The  radar  principle  is  sketched  in  Figure  1.  A  pulse  of  length  ris  transmitted,  is  reflected  at  the  target  and 
is  received  again  at  time  t0  at  the  radar.  From  this  signal  travelling  time  the  range  is  calculated  R0=  ct0/2. 
The  process  is  repeated  at  the  pulse  repetition  interval  (PRI)  T.  The  maximum  unambiguous  range  is 
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therefore  Rmax=  cT  12.  The  ratio  rj=  r/T  is  called  the  duty  factor. 


The  received  signal-to-noise  power  ratio  (SNR)  is 
described  by  the  radar  equation 

R 

Pm 

target  range  (km) 

mean  power  of  transmit  pulse  (W) 

Gt,  Gr 

gain  of  transmit/  receive  antenna 

Oo 

radar  cross  section  of  target  (m2) 

SVR_  1  ps*«_  1  P.G.-^-GA1 

A 

kT0 

wavelength 
=4*10"21  Ws 

(cm) 

(W/Hz) 

(4  xR2f  Pm„,  (4 ttfR'  (kTaFB)-L 

F 

noise  figure 

(dim-less) 

B 

badwidth 

(Hz) 

L 

losses 

(dim-less) 

It  is  the  l/R4  law  that  forces  the  radar  designer  to  increase  transmit  or  receive  energy  as  much  as  possible. 


Fast  time  processing:  The  received  pulse  is  filtered  such  that  the  signal  energy  is  maximally  extracted 
( matched  filtering ,  pulse  compression).  This  is  achieved  by  convolving  the  received  data  samples  zk  with 

L 

the  transmit  wave  form  sk,  k=l..L ,  yk  =  ^Jrzk+r  .  The  range  resolution  after  pulse  compression  is  given 

r= 1 


signal  model  t 


v 


»»»» 


range 
(fast  time) 


by  AT?  >  ct/2,  where  ris  the  effective  pulse  length  after 
compression.  By  suitable  coding  long  transmit  pulses 
with  short  length  after  compression  and  thus  high  range 
resolution  can  be  constructed.  This  requires  a  larger 
bandwidth.  The  ratio  of  the  pulse  length  before  and  after 
compression  is  called  the  compression  ratio  K  which  is 
the  same  as  the  time-bandwidth  product ,  K=  r  befaJ  Tqfter=  B  t  before-  Analogue  waveforms  like  linear 
frequency  modulation  (or  chirp)  are  used  for  pulse  compression,  or  discrete  codes  which  are  switched  at 
certain  subpulses,  e.g.  binary  codes  or  polyphase  codes.  For  radar  the  sidelobes  after  pulse  compression 
are  important  to  avoid  false  targets.  Also,  the  compressed  pulse  must  be  tolerant  against  Doppler 
frequency  shifts  which  a  typical  echo  of  a  moving  target  has. 


Slow  time  processing:  The  receiver  signal  energy  can  be  increased  by  integrating  the  power  from  pulse  to 
pulse.  The  echo  of  a  target  with  certain  radial  velocity  undergoes  a  frequency  shift  fD=  2vR/A  due  to  the 
Doppler  effect.  From  pulse  to  pulse  at  PRI  T  we  observe  thus  a  phase  shift  (f)D=  2  7tfDT.  Maximum  energy  is 
collected  if  this  shift  is  compensated.  The  summation  with  correct  phase  compensation  is  called  coherent 

K 

integration ,  (yj  =  ^  e~jl7lfDkr  yk  .  Of  course,  the  radial  velocity  and  hence  the  Doppler  frequency  fD  is 

k= 1 

unknown  and  must  be  estimated.  The  integration  time  KT  is  called  coherent  processing  interval ,  CPI. 
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I  ^  \  ^  ^ 

Alternatively  one  may  sum  up  only  the  magnitude,  called  incoherent  integration ,  ( |j/|  )  =  34 1  .  The 

'  '  k= 1 

time  of  all  processing  in  a  fixed  look  direction  of  the  radar  (e.g.  with  several  CPIs)  is  called  the  dwell  time. 

2.2  The  phased  array  principle 

The  principle  of  phased  arrays  is  to  generate  a  plane  wavefront  from  a  large  number  of  elementary 
spherical  waves  as  sketched  in  Figure  2.  Some  technical  realisations  of  array  antennas  are  also  shown  in 
this  figure.  The  spherical  waves  are  approximately  realised  by  elementary  antenna  elements  with  near 
omni-directional  characteristic.  The  application  of  suitable  excitations  at  the  elementary  antennas  and  the 
summation  of  all  signals  on  receive  is  called  beamforming. 


Figure  2:  Principle  of  array  antenna  and  some  realisations  of  phased  arrays 


Why  are  phased  array  antennas  of  interest?  The  main  advantage  is  the  nearly  infinitely  fast  switching  of 
the  look  direction  of  the  array.  This  allows  to  illuminate  the  search  space  according  to  some  optimality 
criterion  instead  of  according  to  a  continuous  mechanical  movement.  Recall  that  the  l/R4  law  forces  us  to 
concentrate  transmit  energy.  All  aspects  of  optimising  target  illumination  the  received  energy  are  denoted 
by  the  key  word  energy  management ,  which  is  the  principle  advantage  of  phased  arrays.  Particular 
components  of  energy  management  are 

■  Coherent  integration  can  be  made  nearly  arbitrarily  long.  This  allows  better  clutter  suppression 
(Doppler  discrimination),  target  classification  by  extracting  spectral  features,  and  finally  SAR  and 
ISAR  processing. 

■  Performance  of  different  radar  task  in  time  multiplex,  like  search  and  tracking  of  multiple  targets.  This 
allows  to  use  one  phased  array  radar  as  multi-function  radar. 

■  Optimisation  of  individual  radar  tasks:  optimised  waveforms  for  search,  acquisition  and  track,  high 
precision  measurements  when  required,  variable  beam  shapes,  optimised  algorithms  for  tracking  (the 
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radar  is  “steered”  by  the  tracking  algorithm  and  a-priori  information) 

■  Lower  prime  energy  consumption  (only  for  active  array,  saving  of  about  a  factor  of  2) 

■  High  mean  time  between  failure  (MTBF)  due  to  graceful  degradation  (only  active  array) 

■  If  spatial  samples  over  the  antenna  aperture  are  available:  adaptive  beam  forming  (ABF),  space-time 
adaptive  processing  (STAP),  super-resolution 

2.3  Beamforming 

One  of  the  key  technical  problems  of  phased  arrays  is  the  beamforming  operation.  To  sum  all  signals 
coherently  the  time  delay  of  the  signal  received  at  antenna  element  at  position  r=  (. x,y,z)T  has  to  be 
compensated.  We  denote  the  angle  of  incidence  of  an  incoming  plane  wave  by  the  unit  direction  vector  in 
the  antenna  coordinate  system  u  as  indicated  in  Figure  3  (sometimes  called  the  “direction  cosines”).  The 
plane  in  green  may  represent  the  aperture  of  a  planar  antenna.  The  formulas  are  also  valid  for  3- 
dimensional  arrays.  The  path  length  difference  between  the  element  at  position  r  and  the  origin  is  then  just 

rru=  xu  +  yv  +  -  |r|  •  |u|  -cos^(r,u)  .  For  a  linear  antenna  as  indicated  in  the  left  hand  subplot  of 


=i 


Figure  3  this  is  equal  to  x  sin  9 . 


Figure  3:  Delays  and  direction  cosines 


The  signal  at  element  r  is  written  as  sr  (7,u)  =  b  •  e72^1*  u/c  ,  where /is  the  transmit  frequency  and 

c  the  velocity  of  light.  Correspondingly,  we  form  a  beam  into  a  direction  u0  with  N  antenna  elements  by 
compensating  these  delays 


(1) 


The  superscript  H  denotes  conjugate  transpose.  We  call  a(u0)  the  steering  vector.  For  the  special  case  of  a 
linear  antenna  with  equidistant  elements  at  Xk=  kd/ 1/2  (the  elements  are  separated  by  d/ 1/2)  this  results  in 
the  well  known  function 
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St(u,u0)  =  {t,u)eilKfx^'c  =  Yj->-ellnf'e 


J2*ftjTXk(u-u°] 


k= 1 


k-\ 


N  ,2n  ,  v 

■■b-eJ2aft^e  1  Xl “  “° 


demodulation  k=\ 


j^xk(u-u0)  b  smN7rd(^u-u0) 
d(u  -M0) 


sin^i 


f(u-uQ) 


(2) 


The  function  /(u)=  a(0)//a(u)  is  called  the  antenna  (directivity)  pattern.  From  (2)  one  can  see  that  for  a 
uniform  linear  antenna  this  is  obviously  a  periodic  function,  as  for  any  integer  p  we  have  an  up  such  that 
7tNd  u+27rp=7iNd  up .  This  means  for  element  distances  d>  1  we  have  angles  in  the  visible  region  |  up\<\ 
where  the  main  beam  is  repeated,  called  the  grating  lobes.  Such  grating  lobes  can  be  avoided  if  irregular 
distances  are  chosen.  Array  with  irregular  spacings  are  often  chosen  to  achieve  a  narrow  beamwidth  with 
few  antenna  elements  (e.g.  minimum  redundancy  arrays,  thinned  arrays). 

The  (directivity)  gain  in  a  direction  u0  of  the  antenna  is  defined  as 

GH/(uo)iypfi/(u-u0)i2^u>  (3) 

/  n 

where  Q  denotes  the  visible  region,  i.e.  Q=  {(w,v)  |  u2+v2<  1}.  With  phase-only  beam  steering,  the  gain  is 
approximately  equal  to  the  number  of  elements,  N.  Thus  to  achieve  high  gain  a  large  number  of 
elements  is  needed. 


The  3  dB  beamwidth.  BW  is  defined  as 


,,  ,  BW 

/K±— — ) 


i/2-|/K)f 


For  uniform  linear  antennas  we 


have  BW=  0.887 7A,  or  0BW=  51A4[°],  where  A  denotes  the  diameter  (length)  of  the  antenna  aperture 
measured  in  units  of  the  wavelength  A.  For  planar  arrays  we  have  approximately  BW«1  /A,  or  9Bw& 
60  !A[°\. 


Optimum  beamforming  depends  obviously  on  the  signal  frequency.  The  beamforming  operation  by 
multiplying  the  signal  with  a  suitable  phase  shift  is  therefore  only  correct  for  very  narrowband  signals.  For 
broadband  signals  we  have  to  apply  at  each  antenna  element  the  true  time  delays  (TTD)  rh  =  rf  u0  /  c  .  This 
is  called  broadband  beamforming  indicating  that  a  “phased”  array  is  always  narrowband. 


What  does  narrowband  mean  in  this  context?  First,  note  that  for  planar  phased  arrays  we  have  a 
frequency-angle  ambiguity  exp[y'2^(/  +  A/) rru / c]  =  exp[y'2;r / rT  (u  +  — u) / c] ,  such  that  —  =  — 

u  f 

^Au~" 

and  —  =  — .  From  this  one  can  derive  a  rule  of  thumb  for  the  angle  error  induced  by  a  frequency 

v  / 

deviation  of  Af.  A6[°] «  ^"[0/®]  at  ,  which  is  usually  the  maximum  scan  angle.  One  can  now  set 

an  admissible  beam  pointing  error  at  the  maximum  scan  angle  to  obtain  the  admissible  bandwidth.  For  a 
typical  definition  of  narrowband  one  sets  A6=6Bw  •  This  mean  e.g.  that  for  a  pointing  error  less  than  3°  we 
have  an  admissible  relative  bandwidth  below  3%. 
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A  broadband  source  can  be  written  as  a  sum  (integral)  of  independent  sources  over  all  frequencies  of  the 
band.  For  a  narrowband  phased  array  such  a  source  looks  like  an  extended  source  in  angle 

5/2  5«0/2/0  5v0/2/0 

f  exp\j2zU‘o+f)TTu/c\p(f)nfdf=  f  f  exp\ j27rf0rT  (u0 +u)/ c]  p(u)nadu , 

-B/2  -5«0  /2/0  -flv„/2/0 

where  rif ,  nu  denote  white  noise  processes.  This  is  important  feature  for  suppressing  broadband 
interference.  Note  that  this  frequency-angle  ambiguity  is  only  valid  for  linear  and  planar  arrays,  but  not  for 
true  3-dimensional  arrays.  For  a  3D  array  u  is  always  a  vector  of  length  1,  no  ambiguous  zlu  exists.  For 
too  large  offset  frequencies  the  gain  of  this  type  of  antenna  simply  breaks  down.  This  means  that  3D 
phased  arrays  have  an  inherent  narrowband  filtering  property. 


2.4  Realisations  of  phased  arrays 

The  summation  of  the  possibly  large  number  of  array  channels  can  be  done  either  by  free  space  ( space 
feed  array)  or  a  network  ( constraint  feed  array).  Most  of  the  before  mentioned  advantages  of  a  phased 
array  are  offered  by  the  so  called  active  array ,  where  each  antenna  element  has  its  own  transmitter  and 
receiver,  highly  integrated  in  a  TR  module.  Figure  4  shows  sketches  of  these  types  of  arrays.  The  highly 
integrated  TR  modules  of  an  active  array  are  a  significant  cost  factor. 


Z//Z 

*  *  *  *  emittii 


emitting  antennas 


phase 

shifters 


phase 

shifters 


(c)  active  array 


Figure  4:  Realisations  of  phased  arrays 


A  planar  array  on  a  conducting  plane  has  only  a  limited  field  of  view,  typically  the  sector  of  ±60°.  To 
overcome  this  limitation  several  approaches  have  been  taken.  The  simplest  solution  is  to  rotate  the  array 
mechanically.  This  does  of  course  not  allow  the  full  energy  management  over  360°.  Next,  four  array  faces 
may  be  combined,  but  this  implies  the  fourfold  cost.  Instead  of  four  planar  faces  one  may  arrange  the 
antenna  elements  on  a  curved  surface.  This  is  called  conformal  array.  However,  the  problem  is  the  same  as 
with  the  four  face  array.  The  effective  aperture  uses  only  a  small  part  of  the  antenna  elements  because  the 
other  elements  are  shadowed.  In  addition,  for  a  fully  coherent  integration  in  beamforming  the  different 
polarisations  measured  at  the  differently  oriented  antenna  elements  have  to  be  compensated.  This  means 
full  polarimetric  reception,  i.e.  two  receivers  per  channel,  and  this  means  roughly  a  doubling  of  the  cost. 
Figure  5  shows  the  power  distribution  in  both  polarisations  of  a  spherical  conformal  array. 
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x/a,0  y^o  x/^0  y^0 

Figure  5:  Power  distribution  resulting  from  illumination  of  a  conformal  array  on  a  sphere  by 

plane  wave 


Finally,  the  array  elements  may  be  arranged  in  a  volume.  To  avoid  interaction  with  the  feeding  lines  the 
polarisation  of  the  antenna  elements  has  to  be  chosen  orthogonal  to  these.  This  requires  a  careful  design  of 
the  antenna  elements.  With  this  constraint  and  some  loss  due  to  the  feeding  lines  given,  the  volume  array 
is  the  only  configuration  that  has  full  hemispherical  coverage  with  full  use  of  all  TR  modules.  Figure  6 
sketches  the  experimental  array  called  Crow’s  Nest  antenna  realised  at  FGAN-FHR. 


(a)  Crow’s  Nest  antenna  with  512  elements 
randomly  distributed 


(b)  Concept  of  radar  with  Crow’s  Nest  antenna 


Figure  6:  3D  volume  array 


2.5  Subarrays 

Beamforming  can  be  performed  in  an  analogue  manner  by  a  combiner  network  or  digitally  after  analogue  - 
to-digital  (AD)  conversion  of  the  signals.  For  stability  and  accuracy  we  prefer  digital  beamforming.  For  a 
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high  gain  antenna  with  thousands  of  elements  AD  conversion  of  all  element  outputs  may  be  prohibitive 
with  respect  to  cost.  Cost  effective  solutions  can  be  obtained  by  subarraying.  First  analogue  summation  of 
the  element  outputs  is  performed  after  phase  shifting,  typically  at  RF,  then  the  reduced  number  of  channels 
is  AD  converted  as  indicated  in  Figure  7. 


Mathematically  subarray  forming  can  be  described  by  a  transformation  matrix  T.  Suppose  the  received 
complex  data  at  the  elements  are  described  by  the  vector  zgCw,  let  phase  shifting  for  steering  the  beam 
and  a  possible  amplitude  weighting  be  described  by  a  diagonal  matrix  D,  then  we  can  describe  the 
summation  of  the  data  into  subarrays  by  z  =  T^Dz  ,  where  T  is  of  size  NxL,  if  we  have  L  subarrays.  T 
contains  in  the  z-th  column  ones  for  the  elements  belonging  to  subarray  z,  and  zero  else.  The  phase  shifting 
and  amplitude  weighting  operation  can  be  incorporated  into  the  matrix  T0=  D^T.  The  index  0  indicates  the 
dependence  on  the  chosen  look  direction  u0.  If  the  subarrays  are  disjoint  (non-overlapping),  then  T^T0  is  a 
diagonal  matrix  and  the  receiver  noise  of  the  TR  modules  remains  uncorrelated  at  the  subarray  outputs, 
because  =  E|T0^nn^T0|  =  T|f  Ejnn^jTo  =  cr2!^ T0  is  diagonal.  A  plane  wave  at  the  elements 

transforms  into  a(u)  =  T0//a(u) .  The  elements  of  a(u)  contain  just  the  antenna  patterns  of  the  subarrays 
formed  into  direction  u0. 


An  array  with  subarrays  can  be  viewed  as  a  superarray  having  elements  at  the  centres  of  the  subarrays 
with  antenna  element  patterns  <5.(u) ,  i=l..L.  The  centres  of  the  subarrays  px>h  py>h  can  be  calculated  from 


the  identity  of  the  patterns  5.( u)  =  /J(u)exp 
amplitude  weighting  at  the  elements  Wk 


j2^f(Px.u  +  Py.v)/c 


,  and  one  obtains  for  a  given 


Z  wkSk 


k  eUt 


(4) 


where  we  insert  x  or  y  for  g,  and  C/z  denotes  the  set  of  indices  of  the  z-th  subarray. 

Note  that  the  spacing  of  the  superarray  elements  may  be  greater  than  T/2,  such  that  grating  lobes  may 
occur.  This  is  not  a  problem  if  we  consider  the  beamforming  operation  performed  only  by  the  phase 
shifters,  because  the  full  array  does  not  have  grating  lobes.  The  problem  is  only  apparent  if  beamforming 
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operations  are  performed  at  the  subarray  outputs,  like  steering  a  beam  by  subarray  weighting  or  with 
adaptive  beamforming  at  the  subarray  outputs.  As  we  have  noted  before,  grating  lobes  can  be  avoided  by 
irregular  spacing.  For  a  fully  filled  array  (i.e.  with  elements  at  A/2  spacing)  this  can  only  be  achieved  by 
irregular  subarray  shapes.  This  is  illustrated  in  the  following  figures.  Figure  8  shows  a  generic  array  with 
902  elements  on  a  triangular  grid  at  A/2  spacing  configured  into  32  irregular  subarrays.  We  will  make 
frequent  use  of  this  array  later. 


Figure  8:  2D  generic  array  with  902  elements  grouped  into  32  subarrays 


The  Figure  9  shows  the  resulting  antenna  patterns,  in  subplot  (a)  the  corresponding  sum  beam  pattern  with 
typical  highest  sidelobe  at  -17  dB,  in  subplot  (b)  the  superarray  pattern  due  to  the  subarray  centres 
(without  the  influence  of  the  subarray  patterns),  in  subplot  (c)  the  32  subarray  patterns,  and  the  sum 
patterns  for  some  small  scan  at  subarray  level  in  subplot  (d).  One  can  see  that  grating  lobes  in  the 
superarray  pattern  are  avoided  due  to  careful  irregular  subarray  shaping,  but  scanning  at  subarray  level 
still  leads  to  increased  sidelobes  and  some  loss  in  gain. 
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(a)  Sum  beam  pattern 


u 


(c)  Subarray  patterns 


u 


(b)  Superarray  pattern 


u 


(d)  Scan  at  subarray  level  by  u=  0  (dashed  line),  0.05,  0.1 


Figure  9:  Antenna  patterns  with  generic  array 


3.0  RADAR  DETECTION  AND  PARAMETER  ESTIMATION 

The  first  operation  in  classical  radar  processing  is  the  detection  of  possible  targets.  This  reduces 
considerably  the  data  rate.  Only  for  those  range  cells  where  a  target  has  been  detected  an  accurate 
estimation  of  the  direction,  range,  radial  velocity  etc.  is  performed.  Possibly  a  second  look  is  taken  with  a 
transmit  waveform  that  allows  better  parameter  estimation  (range,  Doppler  etc.).  Such  a  second  look  is 
called  acquisition  or  confirmation.  The  initialisation  of  a  target  track  may  be  performed  after  confirmation 
using  possibly  again  a  different  waveform  optimised  for  the  estimated  target  parameters. 

3.1  Statistical  test  theory 

In  a  statistical  setting  we  consider  the  radar  measurements  (beam  outputs)  zu..zK  as  random  numbers 
consisting  of  an  expectation  value  ^(0),  which  is  the  signal,  plus  noise,  zk  =sk(Q)  +  nk .  0  is  a  general 
parameter  vector.  The  whole  output  data  vector  z  is  characterised  by  a  probability  density  p(z  ;  0),  e.g  by 
the  normal  (Gaussian)  density  N(s(0),  cr2I)  over  complex  vectors.  We  use  the  notation  z  ~  A(0,R)  for  “z 
is  normal  distributed  with  mean  0  and  covariance  matrix  R”.  Detecting  the  presence  of  a  signal  is  a 
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statistical  hypothesis  test  between  two  underlying  distributions,  e.g. 


hypothesis  7i  z  ~  Af(0,  <J2I) 
alternative  A  z  ~  cr2I) 


(5) 


Such  a  decision  is  made  by  calculating  from  the  data  a  test  statistic  T( z)  which  is  compared  with  a 
threshold  7,  and  one  decides  for  H  if  T(z)</j,  and  for  A  otherwise.  This  is  illustrated  in  Figure  10. 


PD  =  Probjr  >rj\A} 
PFl  =  Prob{T 
P f 2  =  Prob  | T  <  rj  \A  } 


detection  probability 
false  alarm  probability 

=  1  -pD 


Figure  10:  Distribution  of  test  statistic  and  detection  and  false  alarm  probabilities 


The  question  is:  how  do  we  get  a  good  test  statistic  ?  One  criterion  for  the  choice  of  the  test  statistic  is  that 
it  should  give  maximum  PD  while  the  false  alarm  probability  is  below  a  given  value  a  which  we  want  to 
control,  and  this  criterion  should  be  fulfilled  uniformly  for  all  parameters.  Such  a  test  is  called  uniformly 
most  powerful  test  (UMP).  It  turns  out  that  the  ratio  of  the  densities  under  hypothesis  and  alternative  with 
the  measured  data  inserted  (which  is  called  the  likelihood)  has  some  optimality  properties.  If  the 
hypothesis  and  alternative  consist  of  a  set  of  densities,  the  maximum  likelihood  has  to  be  taken.  This  is 
called  the  generalised  likelihood  ratio  test  (GLRT). 


max{y?(z;0)} 

Generalised  likelihood  ratio  T (z)  =  9gA  - -  (6) 

max{/?(z;0)} 


The  theory  of  likelihood  ratio  tests  now  gives  us  the  following  results: 

■  If  s(0)  is  a  linear  subspace  (in  particular  if  it  is  only  a  point)  then  the  GLRT  is  UMP.  The  bad  message 
is  that  for  non-linear  alternatives  and  hypotheses  in  general  no  UMP  test  exist.  In  this  case  we  have 
only  the  following  weaker  result. 

■  The  GLRT  is  asymptotically  UMP  for  K—> go  (under  some  general  assumptions). 

The  parameter  that  maximises  the  alternative  is  called  the  maximum  likelihood  (ML)  estimate  of  the 

signal  parameters. 

For  radar  applications  with  zk  being  the  sum  beam  output  for  successive  pulses  we  have  the  following 

results,  see  [2]  Sect.  2,  3: 

■  If  the  system  noise  is  Gaussian  distributed  with  known  cr 2  and  if  the  signal  has  deterministic  but 
unknown,  amplitude,  phase,  Doppler  and  direction,  then  the  GLRT  yields  the  coherently  integrated 
sum  beam  output  as  test  statistic. 

■  If  the  system  noise  is  Gaussian  distributed  with  known  o 2  and  if  the  signal  has  random  phase  and 


RTO-EN-SET -086 


1  -11 


Fundamentals  of  Signal  Processing  for  Phased  Array  Radar 


Doppler,  then  the  GLRT  yields  the  incoherently  integrated  sum  beam  output  as  an  approximation  of 
the  optimum  test  statistic. 

■  The  detection  probability  is  in  general  a  monotone  function  of  the  signal-to-noise  ratio  (SNR). 

These  results  tell  us  when  coherent  and  incoherent  integration  is  optimum,  and  that  for  the  detection 
problem  it  is  sufficient  to  consider  as  optimisation  criterion  the  maximisation  of  the  SNR. 

3.2  Practical  problems  in  radar  detection 

The  measured  noise  power  is  only  a  constant  value  cr 2  if  only  receiver  noise  is  present.  In  reality  we  have 
all  kinds  of  clutter  which  (although  it  is  often  filtered  out  before  detection)  constitute  a  varying  unknown 
(residual)  noise  contribution.  The  noise  power  cr2  used  in  the  threshold  77  is  therefore  often  estimated  from 
neighbouring  range  or  Doppler  cells  of  the  cell  under  test.  The  aim  is  to  maintain  a  constant  false  alarm 
probability  even  in  a  nonhomogeneous  background  interference  situation.  This  is  called  a  CFAR  (constant 
false  alarm  rate)  detector.  The  problem  is  how  to  estimate  the  interference  power,  in  particular  at  the 
clutter  boundaries  (e.g.  from  land  to  see  clutter)  and  for  multiple  targets,  such  that  no  target  is  suppressed. 
The  selection  of  the  cells  for  CFAR  threshold  calculation  and  the  type  of  averaging  is  important  here. 
Mean  and  median  value,  greatest  of,  and  order  statistics  estimation  procedures  have  been  proposed.  All  of 
these  have  advantages  and  disadvantages,  none  of  these  is  globally  optimum. 

A  second  problem  is  impulsive  interference  (e.g.  the  pulses  of  a  neighbouring  radar).  Only  if  such  pulses 
are  received  through  the  main  beam  they  can  represent  targets.  A  detector  to  discriminate  between 
mainlobe  and  sidelobe  echoes  is  therefore  implemented,  called  sidelobe  blanker  (SLB).  This  is  achieved 
by  adding  a  single  omni-directional  element  to  the  antenna,  called  the  guard  channel,  whose  output  is 

■  .2  |  |2 

compared  with  the  sum  channel.  If  \zmain\  <Tj\zguard\  then  it  is  decided  that  the  echoes  come  from  the 
sidelobes  and  the  echo  is  “blanked”.  The  principle  is  illustrated  in  Figure  11. 


Figure  11:  Sidelobe  blanking  with  main  beam  and  guard  channel 


3.3  Maximum  Likelihood  estimation 

In  statistical  estimation  theory  we  consider  the  problem  of  measured  data  z\,..zK  for  which  we  have  a 
parameterised  model  given  by  a  parameterised  distribution  p(zu..zK;  0  ).  We  are  looking  for  a  parameter 
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that  “fits  best”  to  the  assumed  model.  What  does  best  fit  mean?  Any  estimation  0  of  the  desired  parameter 
0  is  a  function  of  the  random  data,  0^(z) ,  i.e.  a  random  quantity.  We  have  two  desirable  properties  that 
should  be  fulfilled,  unbiasedness  and  efficiency. 


The  estimator  0  is  called  unbiased,  if  E 1 0  \  =  0O ,  where  0O  denoted  the  true  parameter  vector. 


For  unbiased  estimators  the  following  bound  is  valid  ( Cramer-Rao  bound ,  CRB): 


E  { A  -  e0 )  A  -  e0  )r }  ^  t  F'1  (0O ) , 


(7) 


where  the  matrix  F  is  called  the  Fisher  information  matrix  and  has  elements 


*U(Oo)  =  J 


(8) 


An  estimator  is  now  called  efficient ,  if  it  has  minimum  variance,  i.e.  if  it  attains  the  CRB, 


For  non-linear  parameter  spaces  unbiased  estimators  rarely  exist.  The  ML  estimator  possesses  this 
property  at  least  approximately.  We  have 

■  For  non-linear  hypotheses  the  ML-estimator  is  asymptotically  unbiased,  for  K-^> go  (sometimes  this 
property  is  also  called  consistency). 

■  The  ML  estimator  is  asymptotically  efficient,  for  go  or  SNR^oo. 

ML  estimators  and  LR  tests  are  often  quite  convenient  to  calculate.  Because  they  have  only  asymptotically 
optimum  properties  one  has  to  prove  by  simulation  that  e.g.  the  bias  is  for  the  considered  sample  size 
really  small  and  that  the  estimator  variance  is  close  to  the  CRB. 

We  remark  that  an  unbiased  estimator  is  not  always  a  good  estimator.  Forcing  an  estimator  to  be  unbiased 
may  increase  the  variance,  so  there  is  in  general  a  trade  off  between  both.  The  total  mean  squared  error  can 
also  be  a  good  criterion.  Other  more  general  optimality  criteria  exist,  like  Bayesian  estimators,  robust 
estimators,  non-parametric  estimators.  For  iterative  processes  like  tracking  unbiasedness  is  in  general  a 
requirement. 

3.4  Monopulse  estimation 

In  this  section  we  will  derive  the  ML  angle  estimator  for  phased  arrays.  It  turns  out  that  this  is  just  a  well 
established  estimator  introduced  earlier  by  ad  hoc  engineering  arguments,  the  monopulse  estimator. 

Suppose  the  model  for  the  data  received  at  the  array  elements  is  z  ~  A(a(u)Z?,  cr2I) .  For  ML  estimation  we 
have  to  maximise  the  complex  normal  density 
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p(z;u,b) 


i*o*Y 


-exp 


-(z  -  (z  -  a(u)6)/  (j2 


(9) 


This  is  equivalent  to  minimising  ||z  -  a(u  )/?||2 .  The  minimum  over  b  is  a  least  squares  problem,  it  ca  be 
written  for  all  u  as  b  =  (adil'  alu))  a(u)"  z  .  Inserting  this  into  the  mean  squared  expression,  using  a;/a  = 
const  and  dropping  constant  terms,  we  are  left  with  the  maximisation  of  the  sum  beam  power 
u  =  argmax|a(u)z|  ,  also  called  the  scan  pattern.  Note  that  this  is  different  from  the  antenna  pattern,  which 

is  the  response  to  source  surrounding  the  antenna  having  a  fixed  look  direction,  while  for  the  scan  pattern 
we  steer  the  beam  electronically  over  the  scene. 

This  scanning  of  the  beam  would  mean  to  transmit  pulses  in  a  fine  grid  of  directions,  which  is  not 
desirable  for  transmit  power  economy.  We  would  like  to  estimate  the  direction  from  a  measurement  which 
is  already  close  to  the  true  direction.  For  detection,  which  is  the  process  before  angle  estimation,  we  will 
scan  the  observation  space  in  a  grid  of  the  beamwidth  BW.  The  3  dB  beamwidth  is  just  the  resolution  limit 
where  the  scan  pattern  over  two  targets  of  equal  power  produces  a  pattern  with  two  distinct  maxima.  In  the 
following  estimation  process  we  then  want  to  estimate  the  target  direction  within  BW  accurately.  The 
estimation  can  be  done  by  approximating  the  scan  pattern  main  beam  by  a  parabola.  We  need  two 
measurements  to  calculate  the  parameters  of  a  parabola,  three  measurements  to  calculate  the  parameters  of 
a  2D  paraboloid  for  azimuth  and  elevation  estimation.  For  monopulse  estimation  one  uses  measurements 
of  the  function  itself  (the  beam  output)  plus  its  derivatives,  called  the  difference  beams,  using  a  Taylor 
expansion  of  the  function.  The  trick  in  monopulse  estimation  is  that  one  applies  the  Taylor  expansion  to 

the  function  G(u)  =  ln  a(u)  z  .  This  results  in  a  numerically  more  stable  linear  extrapolation  formula. 

The  first  order  Taylor  expansion  of  the  derivative  of  G  at  the  position  of  the  maximum  u  is 

G'(u) «  G'( u)  +  G"(u)(u  -  u)  from  which  follows  u  «  u  -  G'(u)/  G"(u) . 


One  can  now  calculate  that  G'(u)  =  2Re<^ 


af(»)z 

a»z. 


5a 

with  the  notation  aA  =  —  for  h=u  or  v.  This  is  just  the 

dh 


ratio  of  the  measured  sum  and  difference  beam  output.  The  elements  of  the  difference  beam  are 
dxi  =dai/du  =  2n jxiai(u) ,  i=  I..N  (if  we  measure  xt  in  units  of  the  wavelength  A)  and  the  analogous 

weighting  for  the  derivative  with  respect  to  elevation.  This  is  a  beam  with  an  additional  amplitude 
weighting  according  to  the  element  position.  The  second  derivative  is  approximated  by  its  expectation 


afa 


value,  which  is  E{G"(u)}  «  2Rej  — >  =  y  (h,t  =  u  or  v).  This  is  a  constant  value  depending  on  the 


a^a 


antenna  element  distribution,  yx  =-N/^xf  .  We  have  thus  obtained  the  classical  monopulse  formulas 

i=i 

([2]  Sect.  3,  [3]) 


u=u0-yxRx(u0 ) 
^  =  vo-/A(uo) 


where  R,  =  Re 


dfz 

a^z 


for  h=x,  y. 


(10) 


All  beams  are  formed  in  a  look  direction  u0,  R  is  called  the  monopulse  ratio. 

For  a  phased  array  three  different  weightings  have  to  be  applied  at  the  elements  for  the  sum  and  the  two 
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difference  beams  which  is  difficult  to  realise.  Therefore  these  different  weightings  are  applied  as  an 
approximation  at  subarray  level.  With  the  digital  subarray  outputs  we  can  form  arbitrarily  many  beams 
with  the  stored  data.  The  smallest  number  of  subarrays  would  the  be  four,  and  the  weighting  with  the 
subarray  centres  then  corresponds  essentially  to  a  weighting  with  ±1  (subtract  left  from  right  aperture  half 
for  azimuth  difference  beam,  lower  from  upper  half  for  elevation  difference  beam).  This  subtraction  of  the 
two  halves  of  is  the  reason  for  the  name  difference  beam  for  the  derivative  beams.  For  the  generic  array 
considered  before  this  is  equivalent  to  replacing  the  subarray  amplitude  weighting  xL  by  sign(xz).  Subarrays 
with  centre  zero  are  not  used  for  difference  beam  forming,  see  Figure  12. 


(a)  Generic  array  as  4  quadrant  monopulse  antenna 


(b)  Azimuth  cuts  through  sum  and  difference  beams 


Figure  12:  Four  quadrant  monopulse  weighting  and  antenna  patterns 


4.0  DETERMINISTIC  BEAM  SHAPING 

In  an  active  phased  array  we  have  in  general  TR  modules  that  can  not  only  control  the  phase,  but  also  the 
amplitude  by  digitally  controlled  attenuators.  This  can  be  used  to  reduce  the  highest  sidelobes,  e.g.  below 
the  -17  dB  level,  which  is  the  highest  sidelobe  of  the  untapered  planar  array,  or  to  reduce  the  sidelobes  in 
certain  angular  sectors  where  interference  is  expected,  or  to  broaden  the  transmit  beam  to  illuminate  a 
larger  sector  to  save  search  time  for  targets  at  close  ranges.  This  is  the  topic  of  this  section. 

Let  us  first  derive  some  general  rules  for  arrays  with  arbitrary  beamforming  weights  weC^  at  the 
elements.  The  antenna  pattern  is 

/(«)  =  |w//a(u)|2  (11) 

which  we  assume  to  be  normalised  to  a  maximum  value  of  1,  i.e.  w//a(0)  =  1.  The  average  sidelobe  level 
then  is 


SL  =  |L  J|/(u)|2  du  =  |T|  Jw^aa^w  du  =  w^Cw 
n  n 

with  C  =  |^|  Ja(u)a//(u)  du .  Note  that  we  have  averaged  also  over  the  main  lobe.  If  the  integration  is 

extended  to  Q=M2,  then  C^I,  and  we  obtain  (under  the  above  normalisation)  a  general  rule  of  thumb  for 
the  sidelobe  level 
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ORGANIZATION 


SL  ~  w/;w  (12) 

In  particular  we  have  for  conventional  beamforming  (w  equal  to  a(,  .V  with  steering  vector  a0)  that 

sl  =  Yn,  \/n2  =  i/n . 

/-^k= 1 


In  Section  3.1  we  have  stated  that  for  optimum  detection  the  SNR  should  be  maximised.  Which 
weightings  maximise  the  SNR?  If  only  white  (receiver)  noise  is  present,  we  have 

SNR  =  E 1 1  w^s  |2  }/E  ||  w^n  |2  j  =  E  jss77 }  w/ E  jnn^ }  w 

=  >02w//aoafw/cr2w/fw  =^>  w  =  a() 

This  says  that  the  conventional  (steering  vector)  weighting  produces  the  maximal  SNR  with  respect  to 
receiver  noise,  or,  in  other  words,  any  other  weighting  producing  lower  sidelobes  has  lower  SNR  against 
receiver  noise. 


4.1  Aperture  tapering 

Now  we  assume  the  complex  beamforming  weights  wi  =  giej27lfVi  u/c ,  i=l.JV.  The  simplest  way  of  pattern 

shaping  is  to  impose  some  bell-shaped  amplitude  weighting  over  the  aperture,  e.g.  g.  =  cosv  {nxi  /  A)  +  a  , 

_  2 

for  suitable  constants  or  g.  =  e  VXi  .  The  foundation  of  these  weightings  is  only  heuristical.  The  Taylor 
weighting  is  optimised  in  the  sense  that  it  leaves  the  conventional  (uniformly  weighted)  pattern 
undistorted  except  for  a  reduction  of  the  first  n  sidelobes  below  a  level  of  x  dB,  [4].  The  Dolph-Chebyshev 
weighting  creates  a  pattern  with  all  sidelobes  equal  at  level  of  y  dB.  Figure  13  shows  examples  of  such 
patterns  for  a  uniform  linear  array  of  40  elements.  The  taper  functions  for  low  sidelobes  were  selected 
such  that  the  3  dB  beamwidth  of  all  patterns  is  equal.  The  conventional  pattern  is  plotted  for  reference 
showing  that  all  tapering  increases  the  beamwidth.  It  depends  on  the  emphasis  on  close-in  and  far-off  low 
sidelobes  which  tapering  is  preferred.  Another  point  of  interest  is  a  small  dynamic  range  of  the  weights, 
because  at  the  array  elements  only  an  attenuation  can  be  applied.  One  can  see  that  the  Taylor  tapering  has 
the  smallest  dynamic  range.  For  planar  arrays  the  efficiency  of  the  tapering  is  slightly  different. 


Figure  13:  Low  sidelobes  by  amplitude  tapering 
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4.2  Optimum  pattern  shaping 

From  a  signal  processing  point  of  view  it  is  the  interference  power  coming  over  the  sidelobes  that  we  want 
to  minimise.  This  can  be  achieved  by  solving  the  following  optimisation  problem,  [5], 

min J|w//a(u)|  p(u)du  subjectto  w^a0=l,or  (13) 

W  Q 

minwCw  s.t.  w/7a0  =  1 ,  with  C  =  Ja(u)a(u)//^(u)<iu. 

As  before,  Q  denotes  the  angular  sector  where  we  want  to  influence  the  pattern,  e.g.  the  whole  visible 
region  u2+v2<  1,  and  p  is  a  weighting  function  which  allows  to  put  different  emphasis  on  the  criterion  in 
different  angular  regions. 

The  solution  of  the  above  optimisation  is 


w  = 


<c-\ 


(14) 


For  the  choice  of p  some  remarks  are  necessary: 

1 .  For  a  global  reduction  of  the  sidelobes,  Q  =  ju  e  R2 1  u2  +  v2  <  1  j ,  one  should  exclude  the  main  beam 

from  the  minimisation  by  setting  p=  0  on  this  set  of  directions.  In  fact,  it  is  recommended  to  exclude  a 
slightly  larger  region  (e.g.  the  null-to-null  width)  to  allow  a  certain  mainbeam  broadening. 

M 

2.  One  may  form  discrete  nulls  in  directions  Ui,...uMby  setting  /?(u)  =  ^£(u  -u*)  .  The  solution  then  is 

k= i 


w  =  — — — - —  with  P  =  I-A(A//A)  lAH  with  A  =  (a(u1),-*-a(uM))  (15) 

a0  Pa0 

To  avoid  insufficient  suppression  due  to  channel  inaccuracies  one  may  also  create  small  extended 
nulls  using  the  matrix  C. 

3.  One  may  reduce  the  sidelobes  in  selected  areas  where  interference  is  expected.  This  is  shown  in 
Figure  14,  where  the  sidelobes  in  the  lower  elevation  space  have  been  reduced  for  the  application  in  an 
airborne  radar  to  reduce  ground  clutter  reflections,  see  [6]. 
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Main  beam  elevation  +60°  Main  beam  elevation  -60° 


-1  -0.5  0  0.5  1  -1  -0.5  0  0.5  1 


u  u 

Figure  14:  Antenna  patterns  with  reduced  sidelobes  at  lower  elevations  for  different  beam 

pointing  directions 

4.3  Phase-only  pattern  shaping 

Sometimes  the  array  may  not  offer  the  possibility  to  apply  amplitude  control,  or  one  does  not  want  to 
apply  any  attenuation  (typically  for  transmit  pattern  shaping).  In  these  cases  pattern  shaping  by  phase-only 
control  is  of  interest. 

Beam  broadening 

To  reduce  the  revisit  time  in  search  mode  one  may  broaden  the  main  beam  on  transmit  and  apply  multiple 
receive  beams  formed  at  subarray  level.  A  simple  broadening  could  be  achieved  by  dividing  the  aperture 
into  four  quadrants  and  applying  in  each  quadrant  a  weighting  for  a  slightly  squinted  beam.  The  sum  of 
these  four  beams  gives  a  pattern  with  broader  beam.  This  would  correspond  to  a  sort  of  triangular  phase 
variation  over  the  aperture.  More  general,  one  can  apply  any  bell  shaped  phase  variation  similar  to  the 
amplitude  variations  given  in  Sect.  4.1.  It  depends  on  the  acceptance  of  ripple  within  the  main  beam  and 
the  height  of  the  first  sidelobes  which  of  these  phase  variations  is  preferred. 
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Phase-only  nulling 

Phase-only  nulling  has  been  a  topic  of  research  since  many  decades.  The  problem  is  that  the  corresponding 
optimisation  problem  (13)  is  not  analytically  solvable.  It  is  even  not  known  if  a  solution  exists  in  general. 
The  parameter  space  (the  phases)  is  a  non-convex  space  which  creates  problems  for  gradient  methods.  For 
well  separated  nulls  iterative  methods  have  been  suggested,  but  these  converge  only  under  the  assumption 
of  small  phase  variations.  It  is  also  known  that  the  creation  of  two  nulls  in  opposite  directions  u  and  -u 
creates  problems.  Some  important  references  here  are  [7],  [8],  [10],  [11]. 

The  convexity  problem  can  be  overcome  by  formulation  of  the  optimisation  problem  and  the  gradient 
algorithm  on  the  phase  manifold.  This  requires  the  use  of  tools  of  calculus  on  Riemannian  manifolds.  This 
approach  and  associated  Newton  iterations  have  been  formulated  in  [9]. 


4.4  Stability  considerations  for  low  sidelobe  antennas 
Influence  of  errors 

Suppose  we  can  apply  instead  of  the  desired  weights  w0j  only  the  weights  w.  =  (|  w0z.  |  +  A g. )  exp  [  jAfy  ] .  In 
addition,  we  assume  that  we  have  errors  in  the  coordinates  of  the  element  positions 
x.  =  x0i  +  Ax.,  yt  =  yQi  +  Ay. .  The  antenna  pattern  then  is 

N 

f(u, v)  =  X  (I  wo,  |  +  Ag, )  exp  \_jln  (  (xi  +  Axi  )u  +  ( v,  +  Ay,  )v)  -  jAfi  ] 

i= 1 


Assume  that  all  errors  are  uniformly  distributed  in  a  given  interval  |Ag.|  <  AG,  |A^.|  <  AO,  |Ax.|  <  AX 

(otherwise  the  elements  are  declared  defective),  then  one  can  calculate  the  expectation  and  variance  over 
the  uniform  distribution.  For  a  linear  array  we  obtain 


E{f(u)}=f0(u) 


sin  AO  sin27ruAX 


AO 


var 


{/(«)} = 


A  G2 


27ruAX 
sin2  AO 


nil 


-  +  1- 


wn, 


(AO)2 


sin2  2nuAX 
(2ttuAX)2 


One  may  now  plot  confidence  interval  patterns  | E[f  (m)}|  ±  var {/ (i/)} .  By  simulation  we  found  that  99% 

of  the  erroneous  patterns  lie  within  this  interval,  so  this  is  a  good  approximation  (for  a  Normal  distribution 
the  confidence  level  would  be  68%).  Figure  15  shows  this  bound  applied  to  the  generic  array  of  Figure  8 
with  a  15  dB  Gauss  taper  at  the  elements  resulting  in  a  peak  sidelobe  level  of -30  dB.  We  have  assumed 
amplitude  errors  of  ±2.5%  of  the  maximum  taper  value,  and  phase  errors  of  ±2.5%  of  2n.  This  is  a  good 
indication  of  the  sensitivity  of  low  sidelobe  antenna  patterns. 
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tapering:  15dB  Gauss,  2.5%  amplitude  errors,  2.5%  phase  errors 


Figure  15:  99%  confidence  intervals  for  sidelobes  with  15dB  Gauss  tapering  and  2.5%  amplitude 

and  phase  errors 


Influence  of  element  failures 

If  say  M  elements  fail  and  are  switched  off  (wz=0),  we  have  the  distorted  directivity  pattern  (11) 

/( u)  =  f(u)-YJliej2*fxIu,cwi. 


The  M  elements  are  a  random  subset  of  all  N  elements.  The  power  density  then  results  in 


/(«)  =|/(u)|  -2Rej/(u)^.=1e 


j  2n  f  xfu/c 


M  il7rf(xi  -xk)Tu/c 


W,Wr 


Now  assume  that  we  have  low  sidelobes,  i.e.  that  |  /  (u)|<cl  in  the  sidelobe  region.  Then  the  following 
approximation  holds  in  the  sidelobe  region 

|/oof 

Because  of  the  assumed  randomness  of  the  defective  elements  the  contributions  for  Nk  will  cancel. 

This  means  the  average  sidelobe  level  is  determined  by  the  power  of  the  defective  elements 

i/wf-sw- 

In  fact,  it  is  very  important  whether  an  element  with  high  or  low  taper  weight  fails.  We  may  further 

II  ||2 

approximate  each  |wz|2  by  the  average  taper  power  w  /  N  and  obtain  the  rough  approximation  for  the 
sidelobe  level 


|/(u)|  «m||w||2  IN . 

For  example,  for  conventional  beamforming  with  wt=  1  IN  failures  of  few  elements  will  have  nearly  no 
effect,  as  the  sidelobe  level  is  already  high.  For  the  40  elements  antenna  of  Figure  13  we  have 
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I  Olog  10(1  /TV  2)=  32  dB  and  sidelobes  only  outside  +0.6  will  be  affected.  For  the  30  dB  Taylor  pattern  we 

have  10-logl0(||w||2  /  TV)  =  31.3  dB  and  nearly  all  sidelobes  may  rise  up  to  this  level. 
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